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Abstract 

From known phase diagram regions of different model Hamiltonians describ- 
ing strongly correlated systems we deduced new domains of the ground state 
phase diagram of the same model by an unitary transformation. Different 
types of extended Hubbard Hamiltonians were used for the starting point 
and the existence of new stable spin-density wave, charge-density wave, fer- 
romagnetic state and a paramagnetic insulator is demonstrated. The used 
procedure itself is dimension independent. 
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I. INTRODUCTION 



Because of their interesting properties, the strongly correlated systems polarized a huge 
effort in last decades given in specially by developments in high critical temperature super- 
conductivity (Degiorgi 1999), heavy- fermion systems (Mathur et al. 1998), metal-insulator 
transition and layered materials (de Boer et al. 1999). The rigorous results obtained in this 
field are related mainly to one dimensional case being deduced by Bethe anzats and exact 
result for D > 1 are extremely rare. Generated by these state of facts, in this paper an 
interesting procedure is presented that allows the construction of new exact ground states 
(GS) in arbitrary dimensions for strongly correlated systems starting from known domains 
of the phase diagram of a given model in the same dimension. The essence of the method 
is simple. Given an examined Hamiltonian (H) together with its ground state ip g valid in a 
restricted region of the phase diagram R, based on an unitary transformation UW = 1 we 
obtain H' = U H W, ip' g = U tp g . The eigenvalues being unmodified by an unitary transfor- 
mation, ip' represents the GS wave function for H' . For a class of analysed Hamiltonians, 
H can be given however with a coupling constant set { gi } and an operatorial terms set 
{ Oi } as H = J2i 9i Oi. Taking U in such a way to transform the components Oi within 
the {Oi} set (i.e. U = O i2 ), the unitary transformation change in fact { & } to { g\ }. 
Based on this observation, from R, the transformed R' region can be deduced, R' being the 
parameter space region where ip' g is valid. 

In this paper we are using an explicitely given U unitary transformation that transforms a 
spin density wave (SDW) into a charge density wave (CDW) wave function, a superconductor 
into a ferromagnet, a phase separation in spin into a phase separation in charge, and vice 
versa. The {Oi} operatorial set is so fixed to preserve the analysed model in the class of 
extended Hubbard models. The concretely analysed cases are connected to the model used by 
de Boer at al. (1995 a, b) and Montorsi and Campbell (1996) describeing superconducting 
properties, and Strack and Vollhardt (1993) and Aligia et al. (1995) studying different 
density wave phases. Our results underline the presence of new ferromagnetic, SDW and 
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insulating GS wave functions in the phase diagram of the enumerated models. 

The remaining part of the paper is conscructed as follows: in Capter II. we describe the 
family of the model Hamiltonians that we analyse and define our unitary transformation U 
used during the paper. Chapter III. contains the study of different concrete Hamiltonians 
and the results obtained. In Capter IV. we present the summary, and Appendices containing 
the mathematical details close the paper. 

II. THE SYSTEM AND METHOD USED 

For the starting point we are introducing below a class of model Hamiltonians that are 
used in the remaining part of the paper for the exemplification of the used procedure. 

A. The studied class of models. 

We are developing here our study in a class of extended Hubbard models that generically 
can be given in the following form 

H = —ti+ U U + X X + V V + J z J z + — J xy + Y Y + R R 

2 

N N 

+ PP + QQ + nY. n i + £ K-t - n n)- (!) 

3=1 3=1 

In this expression a general A A term represent the contribution of the A operator in the 
extended Hubbard Hamiltonian H, taken into account with the coupling constant A. The 
terms taken into consideration, in order, are the following. 
The kinetic energy contribution has the form 

t = ^ ( Cj a Ci a + c] CT Cj a ) , (2) 

<j,i>,<J 

where q ct are canonical Fermi operators wich describe electrons on a d-dimensional lattice, 
< j,i > denoting nearest-neighbor sites. The Hubbard on-site interaction is taken as 

N 1 1 

UU = £/£Kt " «)(n iA - - ), (3) 

3=1 Z Z 
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where U represents the Hubbard coupling and n ia = c\ a Ci a denotes the particle number 
operator for electrons with spin a on site i. The bond-charge interaction is taken as 

* = E ( C )° Ci ° + C Ia C 3<r ) ( n 3-a + "i,-«T ) ■ (4) 
<i,j>,er 

For a density-density type interaction we have taken into account the 

V = £ k - 1)U - 1), (5) 

<j,»> 

contribution that describes a nearest-neighbor Coulomb interaction. 

Heisenberg type spin-spin interactions follows written for spin densities. These interac- 
tions have been taken anisotropic and are represented by 

I = E S j S l (6) 

<j,i> 

Jly = E (SfSr + 5r5+). (7) 
<?',«> 

The spin operators in this expression are standard, namely = 1/2 (n^ — riji ) , Sj~ — 
t c— t 

A pair hopping term follows given by the Hamiltonian term 

^ = E ( c ]t c il c *i c *t + 4t 4| Cjl CjT ) • (8) 

<j,i> 

After this step three-particle interactions are taken into consideration. A correlated hopping 
term testing the opposite spin occupancy along the bond of the hopping is represented by 

R ^ ' ( Qcr Cj a ) Tlj^—u Tli t — a . (9) 

Correlation effects produced by an empty site or a double occupancy are described by 

P = E [Kt - h(*il - h(rn - 1) + (ni T - h(nn - \ ) - I)] . (10) 

<j,i> z z z z 

Finally, a four-particle interaction takes into consideration a double occupancy - double 
occupancy type correlation effect along different bonds 

Q = E ( n iT ~ \) ( n jl ~ \ ) ( n i\ - \ )( n il - \) ■ (11) 



Furthermore, in the starting Hamiltonian given in Eq. ([!]), /x represents the chemical poten- 
tial, and h is the external magnetic field applied in z direction. 

Concerning the notations from the Hamiltonian H we mention some limiting cases that 
emerges during the presentation. When the Heisenberg type interaction Eq.(^) becomes 
isotropic, the notation of the coupling constant becomes W = J z = J xy and the corre- 
sponding Hamiltonian term is denoted by 

W = £ SjS,, (12) 

<j,k> 

where we have Sj Sk = Sj Sji. + 1/2 ( Sj~ + Si S£ ). Furthermore, in case of V — Y , 
the common notation of the VV + YY interactions becomes SS where 



S — $Z S'j S'k — Sj Z S£ + - ( Sj + S k + Sj S k + ) 



<j,k> <j,k> 

w„.,„„ n c'+ _ J J 



(13) 



where Sf = | ( n iT + n n - 1 ) , 5/ = c] T , 5,- = c n c n . 

In the following paragraph we are presenting some basic wave functions that emerge 
during the paper. 

B. The applied wave functions. 

We are presenting in this Section the ground state wave functions used during the paper. 
The superconductor state is considered in the form of an ^-pairing state as follows 

|*,> = (S + )^\0), (14) 
where | ) is the bare vacuum with no fermions present, and 

N 

s+ = £ e< Mi, ( 15 ) 
i=i 

creates the repairs within the sytem. Here N represents the number of lattice sites, and N v 
the number of //-pairs present. 

The ferromagnetic wave-function is taken as 



\*f) = n n 4^1°)' 

ieB jeB' 

where B and B' are arbitrary disjoint sets of lattice points, wich together built up the whole 
lattice. The number of components in B and B' is considered different. 

The density wave functions considered are the following. The spin-density-wave (SDW) 
state is given by 

\*s) = n n 4 T 4io), a?) 

ieA j&A' 

where A and A' represent sublattices of a given lattice. We mention that during the paper, 
in case of density waves we are working on bipartite lattice with odd and even sublattices. 
In case of charge- density wave (CDW) we have 

l*c> = II C ! T 4| 10), (18) 

i&A 

where the A has the meaning of a sublattice. 



C. The unitary transformation used 

The unitary transformation that we analysed during the paper has the form 

N 

u = 11(^4* + ca). (19) 
i=i 

As presented in Appendix l.,U given in Eq. ([l9|) represents indeed an unitary transformation. 
It transforms the presented wave functions (see Appendix 2. for details) as follows 

U | % ) = | *F ) , 

U\* s ) = \* c ). (20) 

Based on Eq.(|19D the Hamiltonian from Eq.(P can be unitary transformed. This procedure 
may be started with the transformation of the creator and annihilator operators (details 
presented in Appendix 3.). 
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t/ci^t = ( _i)iv-i e -iP <Cii< (21) 

Based on Eg. (pi]) the transformation of the particle number operators can be given. As 
presented in Appendix 4., the following results are obtained 

U U ] = m\ , 

u mi u ] = i - mi ■ (22) 

Using Eqs. (^T|,^) the complete Hamiltonian presented in Eq. ([J) can be unitary transformed 
based on U. For example in the case t — X, P = R = 0, and e lPl = 1 for all i, we obtain 
the following final result (see Appendix 5.) 

H' = - it + U U + X X + V V + J z J z + ^ J xy + Y Y 

N N 

+ QQ + flJ2 n j + ~ h £( n iT - n n ) + C. (23) 

3=1 3=1 

In this relation the renormalized coupling constants obtained after the transformation are 
denoted with a tilde superscript their values being 



t = - 


-t, 


Y 


~ 2 ^ xy 


X = 


-x, 


Q 


= Q, 


u = 


-u, 




= h, 


V = 


Jz 
4 ' 


h = 




J z = 


4V, 


C 


= {n - 


Jxy ~ 


: 27. 







(24) 
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III. THE OBTAINED RESULTS 



Based on the results presented in Chapter II. we are analysing different concrete cases 
as follows. 

A. Ferromagnetic states 

De Boer et al (1995 a) in the study of the //-pairing ground state in the extended Hubbard 
model used the following model Hamiltonian 

H l = -tt + UU + XX + VV + J z J z + -^J xv + YY + RR 

2 

N N 

+ PP + QQ + + /i E ( %t - n U ) » (25) 

3=1 3=1 

obtaining a superconducting //-pairing ground state at R = 0, P = 0, in the following 
parameter space region 

V < 0, 
2V = Y , 
t = X , 



- — < max [ V - — + + — + 1^1; 



2 J — i + ^ + 2^+2|t|], (26) 



2; 4 ' 2 

^ 4 

the ground state energy obtained being 



Bw= ^ + ^(V-f + g) + Ar.(2, + ^), (27, 

where iV e represents the number of electrons within the system. 

Based on this result, taking into consideration that from Eq.(^) the unitary transformed 
77-pairing state represents a ferromagnetic state, using Eqs.(^,^) we obtain a ferromagnetic 
ground state in the same model in the parameter space region P = R = and 



S 



Jz < 0, 
z 2 Jxy i 

t = X , 

— < max \ — - V + - \ u \ ; — + V + \Y \ ; 

z 4 z 4 

+ f + y + 2|t|]. (28) 

The ground state energy of this ferromagnetic state is given by 

The position of the obtained ferromagnetic ground state in the T = phase diagram of the 
model relative to the ^-pairing state is illustrated in Fig. 1. 

The second Hamiltonian that we analyse concretely here has been studied by Montorsi 
and Campbell (1996). The authors being interested in the emergence possibilities of re- 
pairing superconductivity, used the Hamiltonian 

N 

H 2 = -ti+Uj2 n iT n ii + X X 

+ 2V S'f S' k z - W W - S ]T S'j S' k + C . (30) 

<j,k> <j,k> 

This Hamiltonian has been studied by Montorsi and Campbell (1996) obtaining repairing 
superconducting ground state for V — in the domain 

U < -A, 

S > 0, (31) 

where A = Z ( 4 | X \ + | | W \ ), Z representing the number of nearest neighbours. 

In Appendix 6. we convert the Hubbard term from Eq.(|3~0|) into the form it has in @. 
The unitary transformation presented in Chapter II then gives a ferromagnetic ground state 
for B.2 in the region 

U < A', (32) 
W > 0, (33) 



where A' = Z ( 4 | X | + | | S" | ), and V = 0, t = X . The position of the ferromagnetic 
domain in the phase diagram of the model is presented in Fig. 2. 

B. Density wave states 

We are analysing now model Hamiltonians giving density wave states. We start the 
presentation of our results with the Hamiltonian used by Strack and Vollhardt (1993) 

H 3 = -tt + UU + X X + V n i° n ™' ■ ( 34 ) 

<j,i>,a,<r' 

The last term from H 3 can be easily transformed into the form used in Eq. (§) (see Appendix 
7.). For the model from Eq . (|3~4"D a paramagnetic insulator has been deduced in the pase 
diagram domain 

t = X , 

|-4t-y>0, (35) 

and a charge density wave region has been obtained in the region (Strack and Vollhardt, 
1993): 

t = X, 

— + At - V < 0. (36) 

Zj 



By the unitary transformation presented above we obtain instead of the phase from Eq.(|35"D 
an other paramagnetic insulator phase present in the region 

t = X, 

-- + 4t - V + > 0, 

Zj 

-tf + 4i + V~ > 0, (37) 

Zj 

and a spin density wave phase instead of the domain Eq. (|36|) in the region 
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u 

- + At + V + > 0, 

Zj 

- - - At + V~ < . (38) 

Zj 

The results from Eqs. (|37| p8|) are valid in a model similar to Eg . fl34|) in which however the 
density-density type coupling has been replaced with V + V + = V + S<j,i>, CT n ja n ia , and 
V-V- = E <j,i>,a n ja n i,-(T j see Appendix 7. The results connected to the Hamiltonian ii/3 
are summarised in Fig. 3. 

The next Hamiltonian that we analysed is connected to an extended Hubbard model 
studied by de Boer et al. (1995 b) 

j „ N 

H A = - tt + U U + X X + V V + J z J z + J xy + Y Y + n Y, rij ■ (39) 

2 3=1 

Based on this Hamiltonian, de Boer et al. (1995 b.) proved the existence of a SDW ground 
state in the following region of the parameter space 

J z > 0, 

Jxy , 



t = X 



7 > max(^ - V; 2 \t\ - ±- V + \Y\ - ±) . (40) 

Using the presented unitary transformation, from Eq . fl40|) we obtain a CDW ground state 
in the following region of the phase diagram of the same model 

V > 0, 

Y = 0, 
t = X, 

_ U > mGX ( - L ]2 \t\ - 2\V\; ^ + \ - V). (41) 
z ~ V 4 16' 11 1 1 ' 2 '2 7 v; 

de Boer and Schadschneider (1995 b.) obtained also a fully polarized ferromagnetic state in 
the following region 



Jz ^ I Jxy | j 

§ > max(2|t| - h + y + il-y _ y + £ ) . ( 42 ) 

Interestingly, from the phase presented in Eq.(^) our unitary transformation gives an insu- 
lating ground state in the domain 
V 

- > 2 Y , 

u . , ol V V J z J z J xy V 

> max( 2 \t H ; — — — ; h — - + — . 43 

z ~ v 1 1 8 16 4 ' 4 2 16 ; K J 

The obtained results related to the phase diagram of the Hamiltonian H4 are summarized 

in Fig. 4. 

Another extended Hubbard Hamiltonian containing in it's ground state phase diagram 
density waves has been analysed by Aligia et al. (1995) 

N 

H 5 = + U n jT n jl + V Y n i° n i<>' + 

j = l <j,i>,aa' 
( c]a C iv + C la C j<? ) [ *AA ( 1 ~ n i-a ) ( 1 ~ n j,-a ) + t B B "^i-a ] • (44) 

<j,i>,<j 

They found a region of parameters in the phase diagram where a CDW phase is the ground 
state of the system in the domain V > V c , where 

V e = , if U < - 2 Z t , 
V c = — + t, 2/ -2Zt < U < 2Zt, 

Vc=^, if U > 2Zt. (45) 

Based on this result, using the presented unitary transformation (see Appendix 8.) we have 
obtained a SDW ground state in the domain V + > V c , V~ < V c , where 

V c = , if - U < - 2 Z t , 
V c = ^— + t if -2Zt < -U < 2Zt, 

V c = -| 2/ -17 > 2Zt. (46) 

In Eq.(p)D we have t = \ t^A I = | £rb| The position of the deduced CDW phase in the 
phase diagram of the model is presented in Fig. 5. 
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IV. SUMMARY 



In this paper, we have presented a procedure that based on an explicitely given unitary 
transformation allows the rigorous deduction of new phase diagram domains starting from 
known ground state solutions of a given model. The method is not model-dependent and can 
be applied for arbitrary dimensions. We have applied this method to find new stable spin- 
density wave, charge- density wave, ferromagnetic and a paramagnetic insulator domains of 
different extended Hubbard type models. 
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FIGURES 



Fig.l. Phase diagram regions in the coupling constant space connected to Hamiltonian 
Hi. a) the superconducting state deduced by De Boer et al. (1995 a.); b) the ferromagnetic 
ground state obtained in this paper. 

Fig. 2. Ground-state phase diagrams for the Hamiltonian H 2 . a) the superconducting 
state of Montorsi and Campbell (1996), b) the ferromagnetic ground state deduced in this 
paper. 

Fig. 3. Phase diagram regions in the coupling constant space connected to Hamiltonian 
H 3 . a) the ground states deduced by Strack and Vollhardt (1993), densely dotted phase is 
CDW state, thinly dotted domain is paramagnetic insulator; b) ground states deduced in 
this paper: the densely dotted phase represents an SDW state and the thinly dotted region 
is an other paramagnetic insulator region. 

Fig. 4. Ground-state phase diagrams for the Hamiltonian H 4 . a) the SDW state obtained 
by de Boer et al. (1995 b.), b) the CDW phase obtained in this paper. 

Fig. 5. Phase diagram regions in the coupling constant space for the Hamiltonian H 5 . a) 
the CDW phase deduced by Aligia et al. (1995), b) the SDW ground state deduced in this 
paper. 
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V. APPENDICES 



In this Chapter dedicated to Appendices we are presenting the mathematical details of 
the calculation. 

A. Appendix 1. 

This appendix prove the unitary nature of the operator U. We must deduce U W = 1, 
where W is the adjoint of U. Indeed, U can be written as 

N 

U = X\ Zi = ziz 2 ... z N , (47) 
i=i 

where zf = ( e~ lPl + ). With this notation we have W = z^ . . . zf, where 
z~l = ( e~ lPl Cii + ). Then we find U W = z\ . . . zn z% . . . zf . Since 

( e iPi c\ { + Cil ) ( e~ lPl Cil + c\{) = 4 Cil + Cil c\ { = 4 Cil + 1 - c\ lCil = 1 , (48) 

we obtain Zi zf = 1 for all i. As a consequence U W = 1, so U is indeed an unitary 
transformation. 

B. Appendix 2 

In this appendix the effect of U on the wave-functions is presented. Let consider first the 
unitary transformation of the superconducting (77-pairing) wave-functions. If we consider a 
single doubly occupied state present in the wave function at site i, the effect of the operator 
Zi on this state is to transform it in a single occuped state with spin up as follows 

Zi\0, Ui, 0, 0) = |0, Ti, 0, 0). (49) 

The effect of Zj on the same state for i 7^ j becomes 

--, 0, '!, .() 0) = e ip i I , . . . , Lj , . . . , Ui , , . . . , ) . (50) 
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From these equations it can be seen that the effect of U on empty sites is to introduce a 
sigle occupancy with spin down on these sites. However, in case of doubly occupied sites 
containing an r/-pair, the effect of U is to transform these states in single occupied states 
with spin up. Based on these relations, for a general wave function containing N v ?7-pairs 
distributed on different sites we obtain 



u\o, o, Un, o, o, UiN v , o, o) = 

T , T , • • • , in , T , • • • , liN v , • • ■ T ) • (51) 



U N e iPj 



Apart from a phase factor, this is a ferromagnetic state, if N v ^ A/2, where N represents 
the number of lattice sites. 

The study of density waves starts from the observation that Zi transforms differentially 
the single occupied states with up and down spin. Indeed, for an up spin single occupied 
state we have 

Zi | • • • , Ti , • • • > = e iPi | . . . , tU , . . . ) , (52) 
but for the spin down counterpart we obtain 

Zi | . . . , U , . . . ) = | — , 0* , — > . (53) 
Using this property, starting from a spin-density wave state, using U we find 

N/2 

U\ |, t, I, !,-••) = II ^ I , Tl , , Tl , • • • ) • (54) 

i=i 

The obtained state is a charge- density state. 

C. Appendix 3 

The transformation under Zi of the creation and annihilation operators is presented in 
this Appendix. Starting with annihilation operators, we have for the spin up case 

Zi Q t = ( e tPl + cn) c n = e lPl q T + c n q T = - e lPl c iT - c n c u , (55) 
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from where 

Zi c iT zf = ( - e lPl q t 4 t - q t Ci A ) ( e~ %Pl Cjj, + ) 

= - c *t c a c *l ~ c *t c *i c l| = - ( 56 ) 

arise. For different indices i ^ j we find z 3 - c^j zj = — q^. Similarly, for the opposite spin 
direction and different indices we obtain Zj z^ = — c^, and for i — j one have 

Zi ai = ( e iPi 4i + en) c u = e iPi c n , (57) 

Zi c tl zf = e iPi ^ c tl ( e~ iPi Cll + cl) = ^ c ti ^ = e iPi ^ (58) 

In the case of creation operators, for up spin, independent on indices we have Zjc\^Zj~ = — 
. For down spin, in case of different indices we obtain Zj c\, Zj~ = — cj, . However, in case 
of i = j the transformation gives 

* c n = ( eiPi 4i + c u ) 4i = c n 4i , ( 59 ) 

Zi c\ z+ = ca cj^e-^cu + cjj = e~ iPi c u Ci[ = e~ iPi Cjj (60) 
^From the presented results the transformation formulas from Eqs.(^) arise. 

D. Appendix 4 

The unitary transformation generated by U of the particle number operators is presented 
below. We have 

U 4 T c lT C/t = U cJ T C/t u ^ = ( _ i c j T ( _ x )iv Qt = c j t = ^ ? (61) 

U cl Cil = U cJi C/t [/ Cii [/t = ( _ ! )iv-i e -iPi c . ; ( _ ! )N -i e P* ^ = 

= ^ cjj_ = 1 - cjj_ Cfj = 1 — n a . (62) 

^From Eq.(|62[) remains unchanged, while is transformed in (1 — n^) under [/. 
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E. Appendix 5 

In this Appendix we are presenting the unitary transformation of different terms from 
the Hamiltonian. We start the presentation with the Hubbard on-site interaction 

U(n n - l/2)(n n - 1/2)17+ = ( n n - 1/2 )( 1 - n n - 1/2) = 

( %T - 1/2 ) ( 1/2 - n n ) . (63) 

It can be seen that the Hubbard interaction term changes the sign under the unitary trans- 
formation used. 

Concerning the hopping term, the following results are obtained 

uftrf = f} 3 

Ufjjrf = -U(c] lCll + 4 lCjl )tf = 

- [ ( - 1 ( - 1 e iP > c n cj, + ( - 1 e~ iPi ( - 1 c n c), ] 

= -e iP(i - i)r f . (64) 

The transformation of the X operator yields 

U Xl C/t = U( ct T Q| + 4 T c n ) ( n n + n tl ) C/t = 

= (cJ T Ci T + c} T c jT )(l - + 1 - rm) = - Xl - 2%. 
U 4 C/t = ( e i(PJ-Pi) Cjl c l + e i(«-Pi) Qi c t i } ( n . T + } 

= - e ^-0 A^ . (65) 

We consider now the J z and V operatorial terms 

U J z C/ f = J C/ ( n iT - riji ) ( 7ii T - n u ) £/ f 
= i ( n jT - 1 + njj. ) ( n iT - 1 + n iX ) = 
_ _ i) = 

[/ y C/t = u ( Wj . _ ! ) ( n . _ ! ) ij\ = u ( njT + Wji _ ! ) ( + ^ _ i ) j/t 
= ( n n - n j4 , ) ( n A - n^) = 4 J z (66) 
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^From Eq. (|66"D the J z and V operators transform into each other. 
The following terms analysed are J xy and Y. We obtain 

U J xy ij = U l -{{S+ Sr + Sr S+) U ] = 
U 2 ( c 1t c n c u c *t + c )i c jt c 1t c ^i) f/t = 

2 ( c lt e ^ c ll e " jPl c *i c *t + c n ^ 4 T e lPl ) = 

1 ( J(Pj-Pi) J J - , - . i pt(Pi-Pi) t r t r .\ _I p i(Pj-Pi)y 

2 ^ e S'T c jl c *i C *T + e c i\ c i[ c 'ji c jV ~ 2 

[/ Fy C/ f = E7 ( c] T Ci| C iT + cJ T Cy; Cft) U ] = 

j(Pi-Pj) t t i i(Pj-Pi) t r ._ r t r ., _ 2( =*( p i- pi ) / .. rfi^l 

e C jT c Ji C i| C *T + e S'J- C ^T C i\ C 4 — ze • \P' ) 

^From Eq . (|67|) , the J xy and F operators transform into each other if the phase factor Pi = 
is considered. Because of this result and Eq.(^), the W and S operators transform into 
each other as well. Indeed, 

TT Q z Q 1 z Tj] Q z Q z 

U Oj U — Oj , 

U S'+ S'r U ] = e i(Pj ~ Pi) Sf Sr , 

U S'f S'+ U ] = e i(Pj - pi) Sr Sf , (68) 

and 

us z 3 stu ] = s'f s'/, 

U S+ Sr U ] = e i{Pj ~ Pi) S'+ S'r , 

U Sj S+ U ] = e iP{i - j) S'f S'f . (69) 

We consider now the transformation of the R operator 

U R)^ = U ( c] T q t + cJ T c jT ) n jl nu = 

( cjt q T + c !t c iT ) ( 1 - nji ) ( 1 - nit) = 

Rji - ( cjt c *t + c !t c iT ) - ( c iT c *t + c !t c jT ) ( n 3i + = 
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U R l Jt lP = U( c] l c ti + c\ t c n ) n n n t] = 
I e Cj| -+- e C4 c^j n, T — 



e i(Pi-PJ) rI. _ ( 70 ) 

We consider now the operator P. This operator is build up from a sum of P^ terms whose 
expression is given by 

1 . 1 . . 1 1 

p ij = ( n n ~ 2 ) ( n 3i ~ 2 n * T + nil ~ 2 ~ 2 > + 

( - ^ ) ( "ti - \ ) ( "it - \ + n Ji - \ ) ■ ( 71 ) 

This operator can be decomposed in the following components 

Hj = ( n ji - 2 ) ( n n ~ \ )( n *i- 2 ) + ( n *T - 2 ) ( n ^ ~ \ )( n n~ 2 ) ' 

^7 = ( n iT - 2 ) ( n Jl ~ ^ ) ( - ^ ) + ( ni T - ^ ) ( n H ~ \ ) ( "jT - ^ ) • ( 72 ) 
The unitary transformation of P + and P~ operators gives the following results 

= (n iT - \ )(-n jir + ^)(-rm + ±) + 

+ (^T " 2)(-^ - 2 )(_n ^ + 2 } = 
C/p.-C/t = ( n . T _ I)(-n j4 . + \ ){n^- 1) + 

+ Kt " ^M-riii + ^)(%T " ^) = -Pij- (73) 

Based on these results one can see, that P£ remains unaltered, while Pj~ changes the sign 
under the unitary transformation used. 

The following operatorial term transformed is Q. We have 

XJQ l3 U^ = U(n n - l -){n n - \ ) ( n A - \ ) ( n n - ±) £/+ = 

( n it - \ )(~ n n + \) ( n n - \) (~ n n + \) = Qa ■ ( 74 ) 

As a consequence, the operator Q transform into itself under the unitary transformation 
generated by U . 
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We consider now the chemical potential // and external field h. Introducing the notations 
fij = rijj + riji and hj = — nji, we obtain 

U fljW = Tlji + 1 — Uji = hj + 1 , 

U hj lr = riji — 1 + riji = fij — 1 . (75) 

We mention that after transforming all terms of the Hamiltonian, driven by the aim to 
remain in the same class of model, we consider the transformed Hamiltonian identical with 
the original. This will be realized if for all lattice sites we have e lPl = 1, and = X^, 
t~l = t\ P+ = p- . Since AT = - X - R and XI = -X, from Eq.(^) we get the 
condition R — . Besides, as = t — 2 X and U = — t we get t = X . We also mention 
that during this paper Pi = has been used. 



F. Appendix 6 

In this Appendix we are presenting the relation that connects the form of the Hub- 
bard interaction analysed in Appendix 4. and the standard form of the on-site interaction 
expressed via a J2i n i,\ n i,i product. We have 

11 11 

12 ( n *T - o ^ ( nil ~ 2^ = ^ ^ Hi] Uii ~ 2 ^ UiX + ^ + 4 ^ 

i i 

= 12 n *t n a ~ \ 12 + 4 12 1 = E ra *T n a ~ 2 ^ + 4" ^ ' ( 76 ) 



il7 



where A e represents the number of electrons within the system and N is the number of 
lattice sites. Denoting A/4 — A e /2 = C, we obtain 

N 1 l 

Ul2n n n n = £/£K T - -)(n n - - ) + UC. (77) 

j=l i 



The result obtained in Eq. ([77|) is used for the study of the Hamiltonian if 2 presented in 
Eq. (|30D . Using Eq. (|77j) , the Hamiltonian H 2 becomes 

N 1 1 

H > = _ tt + C/^(n iT - -)(n n - )+ XX 

j= i * * 

+ 2V J2 S'fS' k z - WW - S ]T S'j S' k + C . (78) 

<j,k> <j,k> 
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G. Appendix 7 

This Appendix presents in detail the transformation of the Hamiltonian H 3 from Eq.( 
in a form more convenient for our study. We have converted the last term of H 3 with the 
help of the following identity: 

Y ( n i° ~ 2 ) ( nia ' ~ \ ) = ^ (( Uja Hi '- a ~ \ Hia ~ \ n i~ a + Z> + 

<j,i>,cr,a' ~~ <j,i>,a ~ 

11 Inn 

Y n ia i n ja - i ]T (n ia + nj _ CT + + n 3 v) + 4 ^ = 

<j,i>,a,cr' <j,i>,a 

Y n ™' n 3° - Y n i° ~ \ Y ( U n + n il + n i\ + n 3i) + 1 = 
<j,i>,<j,a' <j,i>,a <j,i> 

£ n ia . n ja - 2 N e + 1 . (79) 

<j',i>,o-,o-' 



Based on Eq.(|79|) we have 

Y n i°' n i° = Y ( n 3° ~ 2 ) ( ™ i<J ' ~ 2 ^ + COnst - ^ 
Using Eq.([8(]), the Hamiltonian H3 becomes 

H = -ti + UU + X X + VV + const. (81) 



Using the form presented in Eq.(|81|) different terms of the Hamiltonian can be unitary 
transformed. The transformation of the operatorial terms presented in Eq.flSOD gives 

- 2 )( ^ T " 2 } ' 
1 w 1 ^ 

riji )(n i] - ' 



u( 


n n - 








17+ = 




U( 






( n n - 




c/t = 




U( 


n n ~ 




(n iT - 




17+ = 


<|- 


u 


( n ji - 




1 ( n a - 


1 

2 


)C/t = 


<5 



2 
1 

2 



(82) 



Decomposing the contribution presented in Eq. fl80[ ) is spin dependent terms, new notation 
has to be introduced 
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V + = n j° n io > V = U 3° U i-° ( 83 ) 

<j,i>,a <j,i>,a 

The unitary transformation of the operatorial components from Eq. (|83|) gives UV + U^ = V + , 
and UV~ W = - V~. 



H. Appendix 8 



The unitary transformation of the Hamiltonian H§ from Eq.(54) is presented in this 
Appendix. Before starting the unitary transformation we mention that the first two terms 
of ifs may be rewritten based on Eqs.(77, |5D[). For the third term of H$ we are using the 



expression 

■t , . ... J 



12 ( c iT C *T + 4r c it) [ *aa ^1 n jl + t BB ( 1 - thi ) ( 1 - 7lj-i ) ] + 
<i,j> 

12 ( c li c *i + c n c n ) [ ( ( ~ 1 ) tAA ( 1 ~ n n ) ( 1 - ™?t ) + 

<»,i> 

( - 1 ) t BB n A n iT ] . (84) 

The unitary transformation of the new operatorial terms not encountered up to this moment 
are 

U ( cj T q t + c} T CjT ) U ] = ( c] T c iT + cJ T c iT ) 

^ ( 1 - «ii ) ( 1 - n ii) U ] = nn n tt 

U n n U ] = ( 1 - ) ( 1 - n n ) 

u ( c )i c a + 4i c n ) u] = - ( c )i c n + 4i c 3i) 

U ( 1 - n iT ) ( 1 - n iT ) £/ f = ( 1 - n iT ) ( 1 - n n ) 

U ritf %-f Z/t = rijf rij-f (85) 

We mention that during the paper we assumed t^A — ~ 'tBBi an d the following notation 
has been used 

t = I t A A I = I £bb| = I tAA I = I tBB I • (86) 
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